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ABSTRACT 


Using satellite remote sensing data to measure low concentrations of 
suspended solids in lakes and oceans requires careful evaluation of back- 
ground signals from the atmosphere and the v/ater surface. We present here 
typical background corrections for Lake Superior and determine the spect- 
ral distribution of the residual radiance from three major categories of 
turbidity in the lake. The results indicate that for large bodies of 
water, some general information on atmospheric scattering, water clarity, 
and the optical properties of suspended solids allows estimates of con- 
centrations of suspended solids to within ±0.5 mg/Z without using real 
time ground truth data. Under calibrated conditions the threshold detec- 
tion level is 0.3 mg/£ for the fine particulates dispersed throughout the 
laJce and 1 mg/i for the highly light absorbing effluent from rivers. Com- 
parisons of the minimtim reflectance over the open lake areas with reflec- 
tion from the highly absorbing tannin water from rivers, provides a check 
on the clarity of the atmosphere and the excessive background scatter from 
the water surface. 
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METHOD AND RESULT S 

The spectral and the angular distributions of radiance from Lake 
Superior v;ere measured with an optical probe which has a flat response 
from 400 to 900 nrn. The detector had an acceptance cone of 5.65 x 10*'^ 
steradiano and an area of 1 cm^. The spectral distributions of radiance 
were measured using a set of filters spanning the 380 to 1050 nm range. 
The filters had nominal band pass values of 10 nm. Measurements of the 
direct solar intensity were made with the probe and with an auxiliary 
narrow-angle NASA radiometer designed to measure the optical thickness of 
atmosphere. 

The direct solar radiation per unit area normal to the incident 
solar ray is shown in Figure 1. The time of the measurements coincided 
with the LANDSAT 2 overpass. The sun elevation was 57®. The measure- 
ments represent average radiation for June 24, 26/ 27/ and 29, 1979. 

Curve 1 in Figure 1 is based on solar radiation values at the top of the 
atmosphere. The values are based on measurements by Thekaekara (1971) 
and are accepted as standard by NASA (Coulson 1975) . Curve 2 shows the 
radiation at lake level. 

The angular dependence of the radiation from Lake Superior was mea- 
sured at several angles in the incidence plane (the plane containing the 
solar incidence ray and the normal to the water surface). The measure- 
ments v;ere made looking towards the sun (forward angles), away from tlie 
sun (backscattered) , and at right angles to the plane of incidence (side- 
scattered) . Forward scatter is susceptible to glare and is not generally 
used in remote sensing of suspended solids. For light emerging at angles 


near the Zenith, the baukscattered and sidescattered radiances were sim- 
ilar when the sea conditions were calm. The radiances at the Zenith were 
devoid of the excessive specular reflection. However, at angles approach- 
ing 45® from the Zenith, the specular reflection of sky light overwhelmed 
the backscattered signal due to volume reflectance, as shown in Figure 2. 
The intensity and the spectral distribution of the overhead sky light, 
measured at 90® to the solar incidence along a line in tlie plane of inci- 
dence, are sho*n in Figure 3. The overhead sky intensity provides a mea- 
sure of the specular reflection by the water surface. The specular reflec- 
tion must be subtracted from radiance measurements to determine the resid- 
ual radiance from the particulates. Figure 4 shows the spectral distri- 
bution of radiance for calm water containing less tlian 0.2 mg/il of sus- 
pended solids. The radiances are measured near the nadir unless other- 
v/ise specified. The radiance in Figure 4 is broken down into the spec- 
ular component, taken as 6% of the overhead sky intensity, and the remain- 
der of the signal which is attributed to the volume scattering from the 
0.2 mg/il concentration of suspended solids. This remaining signal is 
assumed to be the volume reflectance by the carrying medium. 

In analysis of remote sensing data for turbid water where concentra- 
tions exceed 1 mg/£, the scatter from the carrying medium was diminished 
expor on tially according to the ratio of secchr transparencies. The cor- 
rection term is small, so the use of secchi transparencies is justified. 
Although the secchi transparencies are not precisely measureable, the 
parameter is a convenient and commonly available measure of the transpar- 
ency of lakes and oceans. The decrease in the scatter from the carrying 
medium is important hcicause it implies tliat for highly turbid or highly 
absorbing v/ater, the Ir^ackground is due mainly to atmospheric scattering 
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and specular reflection by the \vater surface. This provides an impor- 
tant checJc on the clarity of tlie atmosphere and the clarity of water in 
the reference area when real time ground truth data are not available. 

For instance, when the minimum signal in the image is high because the 
reference area, assumed to be clear v;ater, actually has turbidity on the 
order of 1 mg/i, the corresponding residual radiance from the kno\im areas 
of highly absorbing tannin water will be unusually low or negative in 
band 4 of the Lcindsat data. In a kno\\m system this v/ill serve as a test 
of the assumption that the minimum reflectance area contains clear water. 
On the other licuid, if the subtraction of the minimum radiance in the 
image yields residual radiances which are high in band 4 over tlie nor- 
mally clear open Icike areas and the areas with highly absorbing water, 
then the excessive background for the day is due to a turbid atmosphere. 

The dependence of volume reflectance on angle is shown in Figure 5. 
The scatter is quite flat over the angles encompassing the range of LriND- 
SAT observation angles and the acceptance angles for the probe, when mea- 
surements were made at or near nadir. 

To obtain residual radiance due to tvirbidity, the measxirements were 
corrected for surface and volume scattering according to Figure 4. As 
previously stated, the volume reflectance from the carrying medium v/as 
reduced exponentially according to the average secchi transparencies. 
Figures 6, 7 and 8 shov; tlie spectral distribution of the residual radi- 
ance for red clay, taconite tailings, and tannin. These 


are the dominant 
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types of suspended solids in Lake Superior. The fine red clay particu- 
lates shown in Figure 6 originate from extensive erosion along the Wis- 
consin shore of Lake Superior. Banks of clay are also present near Onto- 
nagon, Michigan. Erosion of clay is the dominant natural source of inor- 
ganic pcirticulates in the lake. The reflected sunlight peaks at 630 nm 
for the red clay. Red clay concentrations from 1 -15 mg/Z cover large 
sections of the lake, especially in extreme v/estern Lake Superior. 

Taconite tailings give a residual reflectance peak at 560 nm. The 
tailings are discharged from a point source at Silver Bay, Minnesota. 

The discharge, in excess of 6 x 10^ tons/day, causes extensive plumes 
ranging in concentration from 1 to 4 mg/£. Some of the fine tailings 
and the red clay particulates spread throughout the lake and form the 
main background of suspended solids in western Lake Superior. 

Tannin is a term applied to the brown colored river water common 

I around Lake Superior. Tannin is highly organic 35%) . It has low 

transparency, v;ith secchi lengths ranging from .5 - 2 m in comparison 
with 3 - 5 m for western Lake Superior containing 1 mg/Z of red clay 
or tailingc^ The background clear water exceeds 10 m secchi tr'\nsparency . 

i 

1 Figure 9 shows the ^ectral distribution of radiance from rough seas. 

I ■ It was obtained from the difference between radiance measured by looking 

I 

into the wind (corresponding in this case to the sidescatter angles) and 

I the backscattered radiance at right angles to the v;ind. 

The spectral distribution of the residual radiances allow us to 

r 

L identify particulates in the lake (Sydor, Stortz, and Sv/ain 1978). Using 
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the data in Figures 5 and 6 we can estimate the reflectance for each par- 
ticulate type euid concentration. We can also treat mixtures of particu- 
lates, since the reflectances are reasonably linear with tlie concentra- 
tions as shovm for red clay in Figure 6. Hov/ever, to be able to use the 
data in conjunction with the remote sensing information from the satel- 
lites, one needs to be able to predict the radiances at the satellite 
altitude. Furthermore, one needs to account for the seasonal changes in 
solar radiation at lake level and the attenuation of the residual radi- 
ance by the atmosphere. The attenuation by atmosphere can be obtained 
from Figure 1 by taking the ratio of the solar radiation at the lake and 
at the top of the atmosphere. The seasonal variation of the minimum 
background radiance from clear v/ater and the atmosphere is shown by curve 
1 in Figure 10. This variation arises mainly from the effects of solar 
elevation aiid changes of atmospheric attenuation at various solar angles 
(Coulson 1975) . The broken line in Figure 10 approximates the effect of 
solar elevation and is given by 

R = R^(cos3 + asec3)o 

where R is the radiance for clear water and clear atmosphere, 3 is the 
Zenith angle of the sun, a is the scattering coefficient for one atmos- 
phere, and p is the absorption coefficient for one atmosphere, which can 
be evaluated from Figure 1. The limits of fluctuation for background 
radiance are shown by curve 2 in Figure 10. The shape of curve 2 departs 
from the seasonal behavior during the v/inter. This departure can be 
attributed to rescattering. During the summer and fall the albedo over 
the lake is low. Thus correction for atmospheric rescattoring of the 
light reflected from waters adjacent to the observation area can bo 
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ignored. In the winter a high albedo from snow and ice packs surround- 
ing the observed water area accounts for a substantial increase in tlie 
sky light intensity, in the summer and fall curve 2 also corresponds to 
the reflectance for .5 mg/il of suspended solids and minimum atmospheric 
and surface background. In a sense curve 2, for May - November consti- 
tutes the average background radiance for the western Lake Superior. The 
.5 mg/Z concentration is an equilibrium concentration of particulates in 
extreme western Lake Superior. The open lake waters in western Lake 
Superior stay within the .2 - .8 xag/Z limits 97% of the time. Curve 2 is 
the upper limit on the acceptable background signal level for routine 
analysis of remote sensing data in ^-iie absence of tlie real time ground 
truth data. The atmospheric effects are most pronounced in band 4. 

For bands 5 and 6 only the best fit to minimum background values is given 
in Figure 10. Figure 11 shows the direct and diffuse components of solar 
radiation. The curves were obtained by fitting equation (1) to the mini- 
mum values of background radiation for clear atmosphere in bands 4, 5, 
and 6. This determined a, \x, and R^, By scaling to the value of 
solar radiation at the ^op of the atmosphere, we obtain solar radiation 
at lake level. The first term gives the direct component, v;hile the 
second one approximates the diffuse component. The diffuse component is 
not the Scime as the total sky light. 

PROCEDURE FOR REMOTE SENSING DATA ANALYSIS 


When real time sampling data cire available, the background signal 
can be calculated and subtracted from the image to produce residual radi- 
ances for the suspended solids. However in emalysing satellite data 
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without ground truth, we first examine the image to locate the open lake 
area with the minimum signal. If the lowest intensity on the image for 
the open lake areas lies within the 0.2 mg/£ - .5 mg/£ limits shown in 
Figure 10, an assumption of .2 mg/i suspended solids and clear atmos- 
pheric conditions is made. The residual radiance at satellite altitude 
is then obtained by subtracting from the image intensities the minimum 
background intensity given by curve 1 in Figure 10. The resulting data 
can subsequently bo used to identify suspended solids to determine their 
concentrations according to Figures 6, 7, and 8. Hov/ever, if the mini- 
mum signal in the image is substantially above the .2 - .5 mg/£ limits 
shown in Figure 10, a determination must be made of whether the back- 
ground is high because of high overall turbidity of the v/ater or because 
of excessive atrr.ospheric scattering. Visual excimination for extensive 
plumes, thin clouds or haze can be made. Normally, only tapes with clear 
images of the Icike are i^urchased. Hov/ever, a test of excessi.ve back- 
ground signal levels can be made by first subtracting from the image the. 
minimum background according to cufre 1, Figure 10, and then examining 
whether the excess in the background has the spectral character attrib- 
utable to the atmosphere (according to Figure 11), rough seas (Fig. .re 9), 
or suspended solids, (according to Figures 6, 7, and 8). If the spectral 
dependence of the excess in background indicates that the atmosphere is 
turbid, the satellite data may still be usable provided the atmospheric 
turbidity is miform. However, the resulting error for the estimates of 
concentration of suspended solids can be as high as 1 xag/9, and the ident- 
ification of the suspended solids at low concentrations v/ill be hampered 
because of the unpredictable spectral behavior of the turbid atmosphere. 
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The unanibig’.ous identification of the suspended solids relics on 
identifying the ^ectral shape of the residual radiance, and it requires 
a signal at least two or three digital steps above background. This in 
turn requires that the horizontal concentration gradients in suspended 
solids be large enouoh to provide a change of 1 mg/i for red clay and 
tailings and of 3 -ig/'!. for tannin, over the entire area of the lake 
covered by the image. Such .Tadients are readily available for most 
images containing near shore areas and are usually satisfied in the 
images of Great Lakes. However, this condition does Impose some bother- 
some restrictions. In Lake Superior, for example, tannin concentrations 
in the lake must be usually compared v/ith those in the Duluth harbor, or 
a comparable tannin area large enough to be well resolved in the image. 

Although unambiguous identification of particulates from lANDSAT 
data cannot be made for concentrations lower than 1 mg/Z in absence of 
ground truth data, estimates of suspended solids concentration can be 
made to within .5 wg/Z above background. Furthermore, some decisions on 
the identity of particulates can also be made at those concentration 
levels. VJhen no ground truth is available, but the lowest reflectance 
over the lake is within the clearly acceptable background limits, the 
assumption of .2 mg/Z concentratic for the lowest intensity area intro- 
duces a probable error of ± .3 mg/Z in the estimates of concentration of 
red clay and the tailings. The identity of the low concentration of 
suspended solids can also be deduced from spectral dependence based on 
ragnitude of signals in two bands. For v/estern LaJ:e Superior, for 
instance, the background is usually red clay or tailings, so the decision 
on the identity of the low concentration of particulates in the lake can 
be made by exc> mining tlie relative magnitude of the signal in bands 4 and 




ElO 


5. For the tannin water, which usually has higher intensity in band 6 
than band 4, the identity can be based on the relative signal in bands 
6 and 4. Unfortunately, concentrations of tannin lower than 1 mg/i can- 
not be detected in LANDSAT data. At high concentrations, various algo- 
rithms can be devised to automatically extract from the residual radi- 
ances the identity and concentrations of the suspended solids. Figures 
12, 13, and 14 show application of previously published algorithms (Sydor, 
Stortz, and Swain 1978) to the data for June 29, 1979. The identity 
of suspended solids in tlie image agrees closely with the sampling data 
all along the transects of our 1979 cruise of the western Lake Superior. 

COMPARISON OF CALCULATED AND MEASURED RADIANCE AT SATELLITE ALTITUDE 

To exataine how well the results in the previous sections can be used 
to analyse remote sensing data, we calculate the radiances at the sat- 
ellite altitude and compare them with th' satellite data for specific 
sampling sites on Lake Superior. The sampling points are shown in Figure 
15. These points correspond to some of the stations for our 1979 exper- 
imental cruise of Lake Superior using research vessel Crockett . 

The predicted radiance values shotvn in Table 1 are close to the 
observed satellite values. The corrections for atmospheric attenuation 
were taken into consideration according to Figure 1. The satellite data 
were obtained from the Canada Centre for Remote Sensing. The Canadian 
data for LANDSA.T 2 provides a signal resolution level of 1/255. 

Some results discussed in previous sections are evident from Figure 
15 and Table 1. Notice, for instance, that the signal level in band 4 
for the Duluth harbor tannin water is lov/er than the signal in the open 
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lake areas. This indicates excessive background. Thus in absence of 
real time grc'^und truth data, a check on the lov;est background signal 
could have been made to ascertain the clarity of the atmosphere or the 
clarity of the open lake vzater. In this case it would have shown that 
the overall turbidity of the lake \\^s high. The spectral character of 
the residual signal for low intensity open lake areas indicates, accord- 
ing to Table 1, approximately a 1 mg/£ background of tailings. Actually 
the open lake was more turbid than normal, having a background suspended 
solids of level of .7 ng/l of tailings. Thus it is realistic to assume 
that for background conditions falling within the accepted values sho\/n 
in Figure 10, an estimate of the suspended solids could be roade to v/ithin 
± .5 mq/l in cibsence of the ground truth data. Nimbus G data promise to 
provide an even better sensitivity for Lake Superior work. Hov/ever, the 
angular dependence for the volume and surface scattering may be more dif- 
ficult to handle in the Nimbus G data- In anticipation of Nimbus G data, 
it is interesting to compare the spectral shape for the difference in 
volume reflectances observed at a 40® Zenith angle, for 8 mg/£ and 14 mg/£ 
red clay concentrations. This difference is shovm in Figure 16. We 
notice that measurements at Zenith for the 5.6 mg/£ concentration yield a 
similar spectral response. We might thus expect that for Zenith angles 
smaller than ^10® (Figure 2) , v;e might still be a±>le to use remote sensing 
data from the Nimbus Coastal Zone Scanner to identify the particulates 
in plmes, without resorting to elaborate models for cingular scattering 


by small particulates. 
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FIGURE CAPTIONS 


Figure 1 . Direct solar radiation at lake level. The experimental values 
represent 4 day averages. Measurements were raade in .01 pm bands spaced 
at 05 - .08 pm intervals. Curve 1 shows the reference radiation at the 
top of the atmosphere, which is based on standard N.A.S.A. values derived 
from measurements by Thekaekara (1971) . 

Figure 2 . Angular distribution of 0.4 - 0.9 pm light backscattered and 
specularly reflected from clear lake water (.2 - .5 mg/£ suspended solids). 
Measurements were made in the plane of incidence with the instrument 
pointing away from the sun. Specular reflection of sky light becomes 
dominant at Zenith angles exceeding 40° . 

Figure 3 . Diffuse overhead solar radiation measured in the plane of incidence, 
along a line perpendicular to the incident sun ray. The sun was at 23° 

Zenith angle. 

Figure 4 . Curve 1 shows the spectral distribution at nadir of the radiance 
from clearest lake water (less than .2 mg/£ suspended solids). Curve 2 
shows the fraction of this radiance attributable to specular reflection 
of the overhead sky light by the calm water surface. 

Figure 5 . Laboratory determination of the angular distribution of .646 pm 
light scattered by fine red clay particulates and fine mining waste 
particulates. These particulates are the prevalent constituent of suspended 
solids in western Lake Superior. 

Figure 6 . Spectral distribution at nadir of light reflected from red clay 
turbidity in Lake Superior. The residual radiance from particulates is 


linear for concentrations of suspended solids lower than 10 mg/il. 
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Figure 7 . Spectral distribution at nadir of light reflected from mining 
tailings. Tailings often upwell from a broad deposit of discharge slurry 
at the bottom of the lake. 

Figure 8 . Spectral distribution at nadir of light reflected from highly 
organic, opaque river water commonly referred to as tannin. The above 
turbidity contained .5 mg/i of red clay particulates. 

Figure 9 . Spectral distribution at nadir of light scattered from rough seas 
(white caps) . The reflectivity has a flat spectral dependence. Rough 
seas can be easily identified in remote sensing data because of their 
high reflectance at long wavelengths (.8 - 1.1 pm) and unstructured 
geometric patterns. 

Figure 10 . Seasonal behavior of the background radiance due to light 
scattered by the atmosphere, the water surface, and the low concentration 
of particulates in clear water (.2 mg/i suspended solids). The volume 
scatter by the clear water is taken as the intrinsic scatter by the carry- 
ing meditam. Curve 1 represents LANDSAT derived measurements of radiance 
in band 4 (.5 - .6 pm) for the clearest atmospheric conditions. Curve 2 
corresponds to the maximum value of acceptable radiance from clear water 
and uniform atmosphere - defining the clear viewing conditions when 
straightforward analysis for suspended solids may be made in absence of 
real time ground truth data. The broken line in band 4 corresponds to the 
seasonal dependence of radiance as a function of sun elevation. The 
curves for band 5 (.6 - .7 pm) and band 6 (.7 - .8 pm) show average radiance 
for clear viewing conditions. The minimum cutoff radiance in bands 5 and 


6 for LANDSAT 2 and 3 is indicated by the dotted lines. 


Figure 11 . Curves 1, 2, and 3 show the calculated values of direct solar 


radiation at lake level for bands 4, 5, and 6 respectively. Curves A, B, 
and C give the respective total solar radiations in bands 4, 5, and 6, 
including the diffuse component from the overhead sky. The experimental 
points represent measurements of the direct solar radiation plus the 
fraction of forward scattered sunlight included within the acceptance 
angle of the instrument. 

Figure 12 . Distribution of red clay particulate (in excess of .5 mg/£) 
in western Lake Superior. The identifications for suspended solids are based 
on the relative magnitudes of the residual LANDSAT intensity readings in 
bands 4, S, and 6 and the ratios of the residual intensities. 

Figure 13 . Distribution of mining waste tailings in the 1-1.5 mg/£ 
range . 

Figure 14 . Distribution of output from the St. Louis River (lower left 
corner) into western Lake Superior. The output from smaller rivers in 
Wisconsin can also be seen along lower shoreline. 

Figure 15 . Band 4, LANDSAT 2 digital output averaged over 90 pixels. 

Note the low readings in the lower left coiner of the image showing the 
St. Louis River estuary and the Duluth harbor. The overall turbidity 
in the lake is higher than average , thus opaque river water shows lower 


apparent background than the open lake areas. 
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